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Abstract 

We develop a general framework for effective equations of expectation values in 
quantum cosmology and pose for them the quantum Cauchy problem with no-boundary 
and tunneling wavefunctions. Cosmological configuration space is decomposed into 
two sectors that give qualitatively different contributions to the radiation currents in 
effective equations. The field-theoretical sector of inhomogeneous modes is treated 
by the method of Euclidean effective action, while the quantum mechanical sector of 
the spatially homogeneous inflaton is handled by the technique of manifest quantum 
reduction to gauge invariant cosmological perturbations. We apply this framework in 
the model with a big negative non-minimal coupling, which incorporates a recently 
proposed low energy (GUT scale) mechanism of the quantum origin of the inflationary 
Universe and study the effects of the quantum inflaton mode. 



1. Introduction 

This paper is a sequel to the previous work |1J on quantum dynamics of the early Universe 
starting from initial conditions inspired by quantum cosmology. These initial conditions, 
encoded in the no-boundary 0, |3], [|] and tunneling quantum states, can be a source of 
the inflationary scenario at the low (typically GUT) energy scale || [7|, compatible with the 
observational status of inflation theory. In |], [7| the initial conditions were found as a sharp 
peak in the probability distribution of the quantum scalar field whose expectation value 
simulates the effective Hubble constant and drives inflation. The parameters of this peak - 
its location and quantum width - are suppressed relative to the Planckian values by a small 
factor coinciding with the magnitude of the CMBR anisotropy, AT/T ~ 10 -5 , according to 
normalization on COBE || [|. In [10, 11 1]] these results were extended to the open model 



based on the Hawking- Turok instanton 

Although, apriori these features make this model attractive, serious objections usually 
arise regarding the present day status of initial conditions in cosmology. Their issue in 



inflation theory is a subject of strong debate in current literature |13| and numerous meetings 
on structure formation of the Universe. The possibility of quantum cosmological imprint on 
the observational data was called in question either from the viewpoint of the self-reproducing 
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inflation scenario [14], [15[ or from the viewpoints of the eternal inflation [|13fl and anthropic 
principle ||16|| . However, these objections become legitimate only when the conditions of self- 
reproduction and eternal inflation or the dominance of anthropic probabilities are guaranteed. 
In their turn, the realization of these conditions depends on cosmological state. Typically, 
they require the initial Planckian energy scale, a very long duration of the inflationary 
epoch, etc., and as a rule go beyond the scope of reliable perturbative methods. The model 
of m, ||, [7|] avoids all these aspects, because it represents maybe the first example of the 
low-energy phenomenon of the quantum birth of the Universe. Thus, it is likely to be a 
reasonable candidate for the cosmological quantum state deserving further studies. 

The dynamical consequences of the probability peak in the distribution of the inflaton 
were first studied at the level of classical equations || |jj . Then they were reexamined at the 
level of effective equations for quantum expectation values of fields |IJ. However, the analyses 
of [[l] was not complete - the quantum contribution to effective equations included all the 
modes except the main spatially homogeneous mode related to quantum fluctuations of the 
inflaton field itself. This is the main spatially global degree of freedom in cosmology. It is 
responsible for the quantum fluctuations of the homogeneous (minisuperspace) background 
on top of which all the other inhomogeneous modes dynamically evolve. 

As was understood in [|J, the contribution of this mode requires the calculational tech- 
nique different from the technique for inhomogeneous modes. The latter is strongly facil- 
itated by the method of Euclidean effective action which is based on specific properties of 
their quantum state on the quasi-DeSitter background - Euclidean DeSitter invariant vac- 
uum |17] . On the contrary, the quantum state of the inflaton is irrelevant to this vacuum and 
its contribution cannot be obtained by the analytic continuation from the Euclidean effective 
action. Rather, this state is determined by the peak of the inflaton probability distribution 
and can be approximated by the gaussian packet. Direct quantum averaging with respect to 
this packet is required for obtaining the contribution of the inflaton mode to effective equa- 
tions. This step, however, should be preceded by several other calculational steps involving 
the Hamiltonian reduction to the physical sector, setting the quantum Cauchy problem in 
the minisuperspace sector of the system, etc. Thus, the goal of the present paper is to pose 
the quantum Cauchy problem in cosmology for the no-boundary and tunneling quantum 
states, derive the one-loop contribution of the homogeneous mode to effective equations and 
analyze its influence on the inflationary dynamics. 

One of particular questions, to be clarified in this work, is to what an extent this con- 
tribution can change the conclusions of our previous papers [p], 01 • In the spatially closed 
model with the no-boundary and tunneling cosmological states these conclusions were pretty 
stringent. The no-boundary quantum state can give rise only to the eternally long inflation, 
while the finite inflation stage with the conventional exit to the matter-dominated Universe 
exists only for the tunneling stateQ. These conclusions give an undoubtful preference to the 
tunneling state from the viewpoint of the observational cosmology, but, at the fundamental 
level, the tunneling state has an intrinsic problem if one goes beyond the tree-level approxi- 
mationF[ So, there is a hope that the inclusion of quantum fluctuations in the minisuperspace 

Tn open cosmology based on the Hawking- Turok instanton [l2| the initial conditions for finite inflation 
stage can be realized for the no-boundary state jHJ [nj , but the models of (Tl], |l2) with all their merits and 
disadvantages go beyond the scope of this paper. 

2 Point is that the classical Euclidean action enters the exponential of the tunneling wavefunction with 
the "wrong" sign - opposite to that of the loop corrections |L8| . This mismatch invalidates the conventional 
renormalization procedure of absorbing the ultraviolet divergences into the redefinition of classical coupling 
constants |ll| . This inconsistency does not break the results of jl], [| fjj (heavily relying on loop effects) for 
accidental reasons - renormalization ambiguous part of the distribution function in the slow roll approxi- 
mation factors out as an inert overall normalization. However, it shows up beyond this approximation and, 
thus, persists at the fundamental level. 
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sector can render the no-boundary state viable from the viewpoint of the inflationary cos- 
mology. Unfortunately, as we shall see, this hope does not realize - the effect of this mode 
turns out to be too small. This result, however, is model dependent, and the application of 
the general framework of this paper might lead to other conclusions in alternative models 
with well-defined quantum states of the inflaton. 

It should be emphasized, that despite intensive studies on the inflationary dynamics, the 
quantum nature of the homogeneous inflaton mode has not yet been completely understood. 
This mode is peculiar because of its ghost nature - its kinetic term enters the action with the 



wrong sign. This fact was emphasized in [20|, but its cosmological consequences have not 



yet properly been examined. There is a viewpoint that this mode should at all be excluded 
from the quantum perturbations spectrum on the physical grounds that it is always beyond 



the horizon and not directly observable pTfl . This approach can hardly be justified, because 
this mode is in some sense the most fundamental one, for it determines the homogeneous 
background on top of which all observable perturbations are unraveling. The energy of this 
quantum mechanical mode in view of its ghost nature is not positive definite, and it leads 
to peculiar back reaction phenomena which we are going to discuss here. Despite their 
actual small magnitude in our model, in other cases they might essentially contribute to the 
cosmological evolution. As we shall see, quantum fluctuations of the homogeneous inflaton 
have the equation of state p + e = and, thus, it can even be a candidate for the present day 
cosmological constant, instead of quintessence. This is one more motivation for our studies. 

The organization of the paper is as follows. In Sect. 2 we briefly recapitulate the results of 
|], [7], [Tj| for the model with a big negative non-minimal coupling of the inflaton, — £i?<£> 2 /2, 
— £ = |£| 1, having a quasi-quartic potential V((p) = Xip 4 /4 + ... . We formulate the 
mechanism of one-loop radiative corrections due to which the inflaton distribution acquires 
a sharp peak. This peak gives rise to initial conditions for inflation and quantum fluctuations 
contributing to the radiation currents - quantum terms in the effective equations. In Sect. 
3 we reveal the general structure of these currents and their gauge invariance properties. 
In particular, their decomposition into the contributions of two sectors is presented: the 
quantum mechanical sector of the homogeneous inflaton and the field-theoretical sector of 
inhomogeneous modes. The Euclidean effective action method is formulated for 

the calculation of the latter, while the former is supposed to be obtained by direct quan- 
tum averaging with respect to the peak-like wavefunction of the inflaton. In Sect. 4 we 
go over to the generic model with minimally coupled inflaton. For this model we pose the 
Cauchy problem for the classical solution - the tree-level approximation for expectation val- 
ues. We perform quantum reduction of the wavefunction from minisuperspace (of the scale 
factor and inflaton) to the physical subspace, the latter playing the role of the Cauchy sur- 
face in minisuperspace at which quantum initial conditions are posed. For the no-boundary 
and tunneling wavefunctions this subspace is chosen as a domain corresponding to the nu- 
cleation of the Lorentzian quasi-DeSitter spacetime from the Euclidean one. In Sect. 5 a 
similar quantum Cauchy problem is posed for cosmological perturbations. Their physical 



reduction is briefly presented along the lines of pOfl . Perturbations of minisuperspace vari- 
ables - scale factor, lapse and inflaton - are expressed in the Newton gauge in terms of the 
invariant physical variables of |20 |. The latter, in their turn, are found as operators in the 



representation of the physical wavefunction on the Cauchy surface of the above type. In 
Sect. 6 we get back to the non- minimal model of Sect. 2 and discuss its reparametrization 
to the Einstein frame which allows one to transform its quantum Cauchy problem to that of 
the minimal model of Sects. 4-5. We also describe here the calculation of radiation currents 
and derive the expression for the rolling force in the effective inflaton equation. Sects. 7 and 
8 present the resulting effects respectively at the beginning of inflation and at late steady 
inflationary stage. Concluding section contains a brief summary of results and prospective 
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implications of its technique. 



2. Quantum origin of the Universe as a low-energy 
phenomenon 

The starting point of our considerations is the assumption that quantum cosmology 
can give initial conditions for inflation, which in their turn determine main cosmological 
parameters of the observable Universe, including the density parameter Q. The requirement 
of Q > 1 in closed cosmology gives the bound on the e-folding number N - the logarithmic 
expansion coefficient for the scale factor a during the inflation stage with a Hubble constant 
H = a/a, 

N= f tF dtH (2.1) 



o 



(with t = and tp denoting the beginning and the end of inflation epoch). This bound reads 



iV > 60 [O. On the other hand, the value of N is directly related to the initial conditions 



for inflation - initial value of the inflaton, <pi, 

N~- dp^^-. (2.2) 
Jo p 

In the chaotic inflation model the effective Hubble constant H = H(p) is generated by 
the inflaton and, therefore, all the parameters of the inflationary epoch can be found as 
functions of pi. This initial condition belongs to the quantum domain, i.e. it is subject to 
the quantum distribution following from the cosmological wavefunction. If this distribution 
has a sharp probability peak at certain (p — tpi, then this value serves as the initial condition 
for inflation. 

There are two known quantum states that lead in the semiclassical regime to the closed 
inflationary Universe - the no-boundary [Q, ^|, Q and tunneling []5| wavef unctions. They both 
describe quantum nucleation of the Lorentzian quasi-DeSitter spacetime from the Euclidean 
(positive signature) hemisphere - the gravitational instanton responsible for the classically 
forbidden state of the gravitational field. In the tree-level approximation they generate the 
distribution functions 

^.T^-expfT /(¥>)], (2.3) 

where I(p>) is the Euclidean action of this instanton with the inflaton (p. These functions are 
unnormalizable in the high-energy domain <p — > oo and generally devoid of the observation- 
ally justified probability peaks. However, by including quantum loop effects and applying 
the theory to a particular model with strong non-minimal curvature coupling of the inflaton 
one can qualitatively change the situation - generate a sharp probability peak at GUT en- 
ergy scale satisfying the above bound on N |], [7|, |j~8f . The basic formalism underlying this 
result is as follows. 

Beyond the tree level the distribution Pnb,t(<p) is not just the square of the cosmological 



wavefunction (|2.3|) . It becomes the diagonal element of the reduced density matrix obtained 



by tracing out all degrees of freedom but (p. As shown in fll|, |22|, |23], gj] in the one-loop 
approximation it reads 

Pnb,t(v?) ~ exp[=F/(y?) - ri_ loop (^)], (2.4) 
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where the classical action is amended by the Euclidean effective action Pi_i oop (y9) of all 
quantum fields that are integrated out. This action is calculated on the same instanton, 
and its contribution can qualitatively change predictions of the tree-level theory due to the 
logarithmic scaling behavior of r((f). On the instanton of the size 1/H((p) it looks like 
r(ip) ~ Z\nH((f), where Z is the total anomalous scaling of all quantum fields in the 
model. 

The model of || [7j] has the graviton-inflaton sector with a big negative constant — £ = 
|£| ^> 1 of nonminimal curvature coupling, 

S[g^tp] = J *xgW (jj£R(9r) - \^ 2 R(9,u) - \{V V f - ImV - \^ , (2-5) 

and generic GUT sector of Higgs x, vector gauge and spinor fields ip coupled to the 
inflaton via the interaction term 

r A 1 

5 int = J d 4 x g 1/2 ( £ -fx 2 f 2 + £ ^p\ A W + E Uf^ + derivative coupling) . (2.6) 

X A if> 

This model is of a particular interest for a number of reasons. Firstly, from the phe- 
nomenological viewpoint a strong nonminimal coupling allows one to solve the problem of 
exceedingly small A (because here the observable magnitude of CMBR anisotropy AT/T ~ 
10~ 5 is proportional to the ratio \/A/|£| [|25], PB]). Secondly, this coupling is inevitable from 
the viewpoint of renormalization theory. Also, among recent implications, it might be im- 
portant in the theory of an accelerating Universe |[27|| . Finally, for a wide class of GUT- type 
particle physics sectors this model generates a sharp probability peak in Pnb,t(v 9 ) a ^ some 
ip = <pi^, M. This peak belongs to the GUT energy scale - the corresponding effective 

Hubble constant, H(tpi) = yj~\/12\£\(p 1 , is proportional to mpy/~\/\£\ ~ lCT 5 mp. This, in 
its turn, justifies the use of GUT for matter part of the model, because this scale is much 
below the supersymmetry and string theory scales. 

The mechanism of this peak is based on a large value of |£| and the interaction (|2.6|) which 
induces via the Higgs effect large masses of all the particles directly coupled the inflaton. Due 
to this effect the parameter Z (dominated by terms quartic in particle masses) is quadratic 
in |£|, Z — 6|£| 2 A/A, with a universal combination of the coupling constants from (|2.6|) 

^ = ^(E^ + 16£^-i6£4)- (2.7) 

ZA X A V 

Therefore, the probability peak in this model reduces to the extremum of the function 



\np NBtT ((p) ~ const ± - 3— A In — + O — . [2.1 



Here the (^-dependent part of the classical instanton action is taken in the lowest order of 
the slow roll smallness parameter, mp/\^\(p 2 <C 1, renormalization ambiguous parameter /i 
enters only the overall normalization of Pnb,t(v 9 ) an d 



i = (2.9) 
A mi 
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For the no-boundary and tunneling states the peak exists for positive A and respectively 
negative and positive values of 1 + 5, ±(1 + 5) < 0. The parameters of this peak - mean 
values of the inflaton and Hubble constants and the quantum width A, 



37r|l + <5| „, , \/A /27r|l + 5| 

1 %;)=m P ^-J ' (2.10) 



are strongly suppresed by a small ratio \/A/|£| known from the COBE normalization for 
AT /T ~ 10~ 5 [§, Because of small width the distribution function can be approximated 
by the gaussian packet 



Pnb,t(v) - \ /K- ex P 



2A 2 



(2.12) 



The value of the parameter ( |2.9| ) is crucial for the inflationary evolution from this gaussian 



peak. The classical equations of motion in the slow roll approximation, 

+ 3H<p - F{(p) = 0, (2.13) 

show that the inflaton decreases from its initial value, ip ~ F/H < 0, only for 1 + 5 > 0, 
that is only for the tunneling state (minus sign in ( |2.8| )). Only in this case the duration of 
the inflationary epoch is finite with the e-folding number (|2.2|) M 



6ir\£\<p 2 T 48tt 2 

AT~— Jf^- = — — . 2.15 
m 2 P (l + 5) A K ' 

Comparison with N > 60 immediately yields the bound A ~ 5.5 which can be regarded as 
a selection criterion for particle physics models 0. This conclusion remains qualitatively 
true when taking into account the contribution of the inhomogeneous quantum modes to the 
radiation current of the effective equations ||]. This contribution and its dynamical effect 
were obtained in |I[ by the method of the Euclidean effective action, however, the quantum 
fluctuations of the inflaton field itself have not been taken into account. 

For the proponents of the no-boundary quantum states in a long debate on the wave- 
function discord |28|, ^9], [30], |31| this situation might seem unacceptable. According to this 
result the no-boundary proposal does not generate realistic inflationary scenario, while the 
tunneling state does not satisfy important aesthetic criterion - the universal formulation of 
the initial conditions and dynamical aspects in one concept - spacetime covariant path in- 
tegral over geometries^, not to say about intrinsic inconsistency mentioned in Introduction. 
Thus, one of the motivations of considering the quantum mechanical sector of the inflaton 
mode is the hope that it can handle this difficulty In view of the smallness of the quantum 



width ( |2.11| ) the quantum fluctuations Aip ~ A are expected to be negligible, but those of 
their quantum momenta Ap^, ~ 1/A blow up for small A. Therefore, apriori, it is hard to 
predict the overall magnitude of the quantum rolling force and its sign due to A<p(t). In 
what follows we carefully consider this problem. 

3 The Lorentzian path integral for the tunneling state of [fjcfl also requires, in this respect, extension 
beyond minisuperspace level, development of the semiclassical expansion technique, etc. 



6 



3. Effective equations: setting the problem 



Effective equations for expectation values of operators in the quantum state 1^), 

g(x) = (*\g(x)\*), (3.1) 
g(x) = <p(x), x(x),i)(x), Aftix), g^(x), (3.2) 

can be obtained by expanding the Heisenberg equations of motion, SS[g]/Sg(x) = 0, for the 
full quantum field g(x) — g(x) + Ag(x) in powers of quantum disturbances Ag(x) 

S3[g] + r 5 2 S[g] 
8g(x) J 5g(x) Sg(y) 

+ U dvdz s^^lw) AHvmz)+ --- =0 - (3 - 3) 

and averaging them with respect to |^). The linear in Ag(x) term identically drops out of 
this expansion, because (Ag(x)) = 0, and the effective equations acquire a generic form 

Here S[g] is the classical action of the system, and the radiation current J(x) accumulates 
all quantum corrections which begin with the one-loop contribution 

J(x) = \jdydz 8g{x ^ 8g{z) G(z, y) + ... . (3.5) 

The Wightman function of quantum disturbances G(z,y) in a given quantum state |^) 

G(z,y) = (*\Ag(z)Ag(y)\V), (3.6) 
Ag{y)=g{y)-g{y), (3.7) 

can be found by iterations as a loop expansion in powers of %. Because semiclassically 
A(? = Oih 1 ! 2 ) and J{x) = 0(h), it follows from the equations ( |3.3| ) and ( |3.4| ) that the linear 
in A^f term of ( ^.3| ) is at least linear in H. Therefore, in the one-loop approximation the 
quantum perturbation Ag(y) can be identified with the solution of the linearized classical 
equation on the mean-field background 

Its solution can be parametrized by the operator- valued initial conditions. Depending on the 
representation of the initial state they can be either the creation-annihilation operators, 
or operators of initial fields and their conjugated momenta, so that quantum averaging in 
( |3.6| ) becomes straightforward. Continuing this procedure by iterations one can obtain the 
radiation current in any loop order as a complicated but, in principle, calculable functional 
of the mean field. 

Alternatively to (p.5|)-(^6|), the one-loop radiation current can be written as 



J^^ix) 



SS\g] 
Sg(x) 



(3.9) 
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where [...] 2 denotes the quadratic part of the quantity expanded in powers of disturbances 
Ag that solve linearized classical equations, and (...) implies the quantum averaging with 
respect to |^). 

For the cosmological system this simple perturbative scheme should, however, be amended 
by two important aspects. One of them reflects the local gauge (general coordinate) invari- 
ance of the problem and the other deals with disentangling the collective variables. The 
latter describes the most important (minisuperspace) cosmological degrees of freedom hav- 
ing non-trivial expectation values. In the next two sections we consider these two aspects of 
the problem. 



3.1. Gauge properties of the radiation current 

In view of local general coordinate and other gauge invariances, fields and their pertur- 
bations contain purely gauge variables that should be gauged away. Thus, the physical sector 
should be disentangled from the full configuration space of the system and the physical state 
should be prescribed on this physical sector. This is the general scheme of the reduced phase 
space quantization |24], |32], |33|. For describing this scheme in application to perturbative 
radiation currents we simplify the formalism by using condensed notations for the full set of 



fields (|3.1|), g a = g(x), in which the condensed index a includes both discrete spin labels and 
spacetime coordinates x. Contraction of these indices implies also the spacetime integra- 
tion. In these notations, the invariance of the action S[g] under local gauge transformations, 
g a — > g a + R a with infinitesimal gauge parameters / M (the condensed index \i bearing 
together with discrete tensor labels spacetime arguments of the local function / M = f(x)) 
reads 

Kir- = Q - ( 3 - 10 ) 

Here, is a generator of the gauge transformation - the quasilocal two-point kernel with 
respect to spacetime coordinates associated with condensed labels a and /i. 

A gauge breaking procedure - a part of the physical reduction - can be enforced by adding 
to the classical action the gauge-breaking term and, in the quantum domain - for Heisenberg 
equations - by including the action of Faddeev-Popov ghosts. Then the derivation of effective 
equations repeats the perturbative scheme of the above type with gauge-breaking and ghost 
terms included into the full action. However, for the purpose of physical reduction gauge 
conditions should be unitary, that is transforming under gauge transformations locally in 
time (not involving time derivatives of the gauge parameters In such a gauge the ghosts 
are not propagating, and the ghost action can be omitted from the total quantum system. As 
a result, the one- loop effective equations again take the form ( |3.4|) with the same radiation 
current 

1 r) 3c > 

The only modification due to local invariance is that the linearized equations of motion fl3.8Q 
for quantum perturbations Ag a are supplied with the linear gauge conditions 

Ag b = 0, (3.12) 



5g a 5g b 

X »Ag a = 0. (3.13) 
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The functional matrix (two-point kernel) of the linear gauge condition \a ls assumed to form 
the Faddeev-Popov operator 

QZ = xZK, (3-14) 

which is ultralocal in time, Q% ~ i„) (the property of the unitary gauge), and invertible. 
In view of this ultralocality the inverse of Q^ ll/ , does not require imposing any boundary 
conditions in the past or future of the time variable. 

The system of equations ( |3.12| )-( ]3.13| ) for quantum perturbations Ag b has a number 



of peculiarities. First, the linearized gauge condition ( |3.13| ) guarantees that the gauge- 
breaking term (usually quadratic in gauge conditions) does not contribute to the radiation 
current ( |3.11| ). Second, it fixes uniquely the solution for Ag a under given initial conditions. 
In the absence of gauge conditions, eq. ( p. 12 ) would have the ambiguity in the solution, 



Ag a — > Ag a + with arbitrary / M because of a simple corollary of ( 3.10Q 



^ * (3,5, 



5g a 5g b M 5g a 5g b 

Here the right hand side vanishes on the classical solution, 5S/5g a = 0. So, the gauge 
generators are zero vectors of the Hessian of the action on this background, which implies 
the gauge invariance of the linearized solution of the above type. However, the auxiliary 
condition ( p.!3| ) gauges this invariance away and guarantees the uniqueness of the solution 
for Ag a . 

The parametrization of the general solution of eqs. ( 3.12 )- (|3.13|) in terms of the sym- 



plectic (phase space) initial conditions is equivalent to the Hamiltonian reduction of this 
linearized system to the physical sector. This reduction should be done in the canonical 



formalism. The unitarity of gauge conditions ( 3.13|) guarantees that they can be rewritten 
in terms of the phase space variables - configuration coordinates and conjugated momenta 
- contained in the set of Ag a and dAg a /dt (the rest of the variables in Ag a are the La- 
grange multipliers). Solving these canonical gauge conditions together with the first class 
constraints - the nondynamical subset of eq.( [3.12] ) - one finds all the perturbations Ag a 
as functions of the physical variables Ag p h ys which in their turn become functions of initial 
physical coordinates and momenta {q ,p ) 

Ag a = Ag a (Ag phys ), (3.16) 
Ag phys = Ag phys (q ,p ). (3.17) 

If the initial quantum state is known on the physical sector in the representation of these 
quantum variables, |^) = &(q ), p = d/idq , then the calculation of averages in ( |3.11Q 
becomes straightforward. 

In what follows we use this calculational strategy. The no-boundary and tunneling wave- 
functions as solutions of the Wheeler-DeWitt equations on superspace will, first, be reduced 
to the physical sector. This quantum reduction in the one-loop approximation will be per- 
formed by the technique of |24|, |32], [33|[]. Simultaneously, the physical reduction will be per- 



formed for the cosmological background and its perturbations, the both being parametrized 



4 The fact that the cosmological states are known as solutions of quantum Dirac constraints on superspace, 
and the fact that their quantum reduction to the physical sector is readily available from p4| , |3^ , [33| , explains 
why we work within unitary gauge hxing procedure. Relativistic gauges with propagating Faddeev-Popov 
ghosts would require a quantum state on extended Hilbert space with indefinite metric, satisfying the zero 
BRST-charge equation (see review of this problem in pi|). Lifting the Dirac wavefunctions of the no- 
boundary and tunneling states to this space, to the best of our knowledge, has not been done and goes 
beyond the scope of this paper. 
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in terms of initial data encoded in the (reduced) physical wavefunction. This makes further 
calculation of radiation currents straightforward. 

We finish this section with gauge invariance of the radiation current - the corollary of the 
Noether identity for the classical action ( p.lOj ). The latter, after two consequitive functional 
differentiations, yields another identity 

jja S 3 S _ M% 5 2 S 8J£ 5 2 S 5 2 R; 5S 
» 5g a 5g b 5g c 5g b 5g a 5g c 5g c 5g a 5g b 5g b 5g c 5g a ' ^ ' ' 



Using it, one shows on account of the linearized equations ( 3.12|) that the radiation current 
satisfies the relation 



^ji-ioop = _f3_ M { A ^Ar>. (3.19) 
^ a Sg b 5g c 5g a N ' 



Here the right hand side vanishes on the classical background, 5S/5g a = 0, and, moreover, on 
an arbitrary mean field background when the generator is linear in the field, 8 2 R a ^/ 5g b 5g c = 0. 
But this is a well known property of the generators of general coordinate transformations 
that form the closed algebra of spacetime diffeomorphisms. Thus, the one-loop radiation 
currents also satisfy the Noether identity 

R a ji-ioop = Q (3 20 ) 

This property will be very important in what follows. It implies that the radiation currents 
are linearly dependent, which reduces the number of effective equations to be solved in cos- 
mological applications. Moreover, this identity reflects the gauge invariance of effective equa- 
tions themselves. In particular, for purely gravitational system with g a (x) = g^ v (x), when 



-loop i 



X 



the radiation current coincides with the expectation value of the stress tensor 1 
(T Ml/ (x)), this property signifies the covariant conservation law, V At (T Aliy (x)) = 0. 

3.2. Two configuration space sectors of the model 

In closed cosmological model, the total metric and inflaton scalar field are usually de- 
composed into the spatially homogeneous background and inhomogeneous perturbations 

ds 2 = -N 2 (t)dt 2 + a 2 (t)- fij dx i dx j + h^(x)dx^dx v , (3.21) 
ip{x)=ip(t)+8ip(x), x = (t,x), (3.22) 

where a(t) is the scale factor, N(t) is the lapse function and 7^ is the spatial metric of the 
3-sphere of unit radius. Therefore, the full set of fields consists of the minisuperspace sector 
of spatially homogeneous variables Q(t) and inhomogeneous fields f(x) essentially depending 
on spatial coordinates x i= x 

s(z) = Q(t),/(t,x), (3.23) 
Q(t)=a(t),<p(t),N(t), (3.24) 
f(x) = 6<p(t, x), h^(t, x), x(t, x), ip(t, x), A^(t, x), ... (3.25) 

From the structure of the initial quantum state, that will be discussed later, it follows that 
only minisuperspace variables have nonvanishing expectation values 

(Q(t))^0, (/(a;)) = 0. (3.26) 
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Therefore, the full set of effective equations reduces to the following three equations in the 
minisuperspace sector 

SS[Q] + T(t)-Q 

mj + jQit) ~ °' 

Jq = Jn, Ja, Jip, (3.27) 

their quantum radiation currents Jq(£) containing the contribution of quantum fluctuations 
of minisuperspace modes themselves and those of spatially inhomogeneous fields. 

The set of these equations is, however, redundant in view of the Noether identities for 
both classical ( 3.10 ) and quantum ( |3.2U| ) parts. In the minisuperspace sector the general 



coordinate transformations reduce to reparametrizations of time. Infinitesimal transforma- 
tions of minisuperspace variables Q(t) = (N(t), a(t), (pit)), in the notations of Sect. 3.1, 
read as: /** = f(t), = V Q / = (d(Nf)/dt, df, <pf), and the identity ( pT20| ) takes the 

form 

(pJ^ + dJ a - NJ N = 0. (3.28) 

The currents and J v have direct physical interpretation in terms of the quantum energy 
density e and pressure p 

J N = -a 3 ^fye, J a = 3a 2 yfyp, (3.29) 

so that eq. ( 3.28|) in the cosmic time, N — 1, takes the form 



e + 3-(e + p) + ^^ip = 0. (3.30) 



a 



As we shall see later, the third current J v can be interpreted in terms of the quantum 
rolling force driving the evolution of the inflaton field. Therefore, this equation measures 
the balance of the conservation for the quantum stress tensor vs the work of this force. In 
the slow roll regime, ip ~ 0, it reduces to the conventional covariant conservation law on the 
Robertson- Walker background. 

The Hamiltonian reduction to the physical sector, discussed above, 

Ag(x) = (AQ(t), f{x)) -> Ag phys {x) = (AQ phys (t), f phys {x)), (3.31) 
/ P hys(a ; ) = (h TT (x), matter fields), (3.32) 

leaves us with the set of physical variables Ag p ^ ys (x) which also splits into minisuperspace 
and field-theoretical subsets. Here AQ phys (t) is a single field variable that originates from 
the minisuperspace sector of metric and inflaton perturbations^ The rest, f p hy S {%), rep- 
resent the transverse traceless modes of the gravitational wave ^ TT (a;) and other physical 
nongravitational degrees of freedom. The nature of AQphys(^) depends on the gauge used 
for disentangling the physical sector. A particular gauge fixing procedure widely used in 



the theory of cosmological perturbations |2Q, 34] picks up a special gauge invariant variable 



AQphys(t) = q(t) that will be discussed in much detail in Sect. 5. At the nonlinear level (be- 
yond perturbation theory on a classical background) another choice of Q p h ys (£) is possible 
by simply identifying it with the spatially homogeneous inflaton field <p(t), Q p hy S (t) = <p(t). 
In both cases, particular expressions for the original minisuperspace variables Q(t) and their 

5 Counting the number of physical degrees of freedom is usual: in 3-dimensional configuration space of 
(N, a, if) subject to one first class constraint the number of physical degrees of freedom equals 3 — 2x1 = 1. 
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perturbations AQ(t) in terms of Q p h y s(t) and AQ p h ys (£) depend on the choice of gauge con- 
ditions. Two types of these gauge conditions will be considered in Sects. 4 and 5. 

Splitting the whole configuration space into minisuperspace and inhomogeneous sectors 
( |3.23|) reflects the choice of the collective variables. Moreover, it reflects distinctly different 
nature of quantum states for the modes belonging to these two sectors. This results in 
different calculational strategies for the corresponding quantum averages. To see it, note 
that on the space of physical variables ((p, / p h ys (x)), with / p h ys (x) treated perturbatively, the 
initial no-boundary and tunneling wavefunctions read (see Sect. 4) 



* 1 - loop (¥,,/ phys ) = P(^)exp 



^)/ P 2 hys + 0(/ p 3 hy J 



(3.33) 



Here -P(^) is the loop prefactor and the tree-level exponential contains the Euclidean action 
expanded up to a quadratic term in inhomogeneous modes, Q({p)fp hys . 

From (|3.33j ) it follows that the one-loop quantum correlators between the minisuperspace 
modes and inhomogeneous fields vanishP], therefore they contribute additively to the total 
one-loop radiation current 



J =Ji + J f 



j'(t) 



dt' dt" 



dx dy 



S 3 S[Q] 



6Q(t) 5Q(f) 6Q(t> 
5 3 S[Q + f] 



5Q(t) 5 f{x) 5f(y) 



(AQ(t')AQ(t")), 

(mm)- 



(3.34) 
(3.35) 

(3.36) 



/=o 



Similarly to ( |3.9| ) these contributions represent the quantum averages of the quadratic terms 
of the expansion of SS/5Q correspondingly in AQ(t) and f(x). However, the calculation 
methods for J q and J* are very different, and the difference can be attributed to qualitatively 
different quantum states of the modes if and / p h ys - Let us begin with the radiation current 
J* which can be obtained by the effective action method 

The matrix of quantum dispersions fl(<p) in the gaussian state of inhomogeneous modes 
is such that this state turns out to be the Euclidean quasi-DeSitter invariant vacuum ]17| |35 | 



vac) DS = C%)exp 



The corresponding quantum averages 



DS 



(vac|/ phys (x)/ phys (y)|vac) D5 = G DS (x,y) 



(3.37) 



(3.38) 



are given by DeSitter invariant Green's functions which can be obtained by analytic contin- 
uation from the unique Green's function on the Euclidean section of the DeSitter spacetime 
- the inverse of the Hessian of the Euclidean action 



G DS (x,y) = G E (x E ,yE) 
5 2 I[Q + /] 



/ 



dyi 



Sf(x E ) Sf(y E 
I[Q + f} = -iS[Q + /] 



G E (y E , z E ) = S(x E , z e ) 



(3.39) 
(3.40) 
(3.41) 



6 Note that due to the gaussian nature of the state / ~ h 1 ^ 2 , so that the terms contributing to (/A<p)- 
correlators, 0(f 3 ) ~ h 3 ^ 2 , go beyond the one-loop approximation. 
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Here xe denotes the coordinates on the Euclidean DeSitter manifold related by the analytic 
continuation to the Lorentzian spacetime coordinates, x E = ir/2H + ix°, x# = x, with H - 
the Hubble constant or the inverse radius of the Euclidean 4-sphere[| and x° = t - the cosmic 
time in the unperturbed DeSitter metric ( |3.21|) corresponding to N = 1, a(t) = cosh(Ht) / H . 

In view of this relation and in view of a similar analytic continuation rule between the 
Lorentzian, 5 3 S/5Q5f 2 , and Euclidean, 6 3 I/5Q5f 2 , 3-vertices, one finds the Euclidean ef- 
fective action algorithm for the radiation current of inhomogeneous quantum modes 



x 



5Q(x*) 



A-ioo P [Q] = - Trln 



2 5f(x E )6f(y 



E 



(3.42) 
(3.43) 



/= 



Note that this is exactly the functional that yields the one-loop contribution to the dis- 
tribution function (|2.4| ), when evaluated at the classical solution for the minisuperspace 
background Q{x 4 ) = Q(x 4 , tp) parametrically depending on ip, 

A-ioopCv) = JV loop [Q(x 4 ,^)]. (3.44) 

This algorithm was used in the previous paper |I[ for the calculation one-loop radia- 
tion currents of /-modes. In i^i_i 00 p[Q] was obtained by the local Schwinger-DeWitt 
expansion - the expansion in spacetime derivatives of the background fields, which in the 
cosmological context corresponds to the slow-roll expansion^]. Within this expansion the 
one-loop action is represented as a spacetime integral of the effective Lagrangian expanded 
in powers of curvatures, matter field strengths and their covariant derivatives. Therefore, the 
analytic continuation rule ( |3.42| ) is trivial - the current J* (x°) can be given by the functional 
variation of the local Lorentzian one-loop action 

JV) = (3 45) 

Si-tooplQOc )] =iA-ioop[Q(z 4 )]| ++++ ^_ +++ . (3.46) 

Therefore, within the local expansion the J? part of the radiation current can be absorbed 
in the functional variation of the total Lorentzian effective action S e g[Q], and the effective 
equations acquire the final form 

^gW^O, (3.47) 
S eB [Q) = S[Q] + S 1 -i oap [Q]. (3.48) 



7 The Lorentzian Green's function Gds{x, y) = iG^g(x, y) is the positive frequency Wightman function - 
the solution of the homogeneous linearized equation of motion. On the contrary, Ge{xe, Ve), as an inverse 
of the Hessian, solves the inhomogeneous equation. However, the Wightman function can be obtained from 
Ge{xeiVe) by taking the boundary value of its analytic continuation on a proper shore of the cut in the 
complex plane of \ a(x, y )] 1 I 2 - the geodetic distance between the points x and y (TtJ . 

8 The algorithm ( 3.42| )- ( ]3T43 ) looks as a generalization of the analytic continuation method for the effective 



equations in asymptotically flat spacetime | |36| . Strictly speaking, this algorithm as derived above holds 
only for exact DeSitter background, while the method of p6| was proven for arbitrary asymptotically flat 
backgrounds that are perturbatively related to flat spacetime. For the deviations f rom t he DeS itter geometry 



(measured by the magnitude of the slow-roll smallness parameter) the relations ( |3.42| )-( |3.43[ ) hold for local 
terms of the effective action, and might not be true for an essentially nonlocal part that cannot be expanded 
in powers of derivatives. But for the slow-roll inflation regime, which we use throughout the paper, such an 
expansion - the local Schwinger-DeWitt series - is definitely applicable, which justifies the effective action 
method. 
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Here S'cfrfQ] can be obtained from the classical action S[Q] by adding loop corrections to 
the classical coefficient functions in the curvature and gradient expansion of the Lagrangian 
f\. These corrections, in their turn, can be calculated in the Euclidean spacetime by the local 
Schwinger-DeWitt technique [j37fl . For our model these corrections has been obtained in [0. 
Thus, this essential simplification in treating the J^-part of the radiation current exists due 
to the two important aspects of the problem - DeSitter-invariant vacuum of /-modes and the 
slow roll approximation. In contrast to this, to the best of our knowledge, no simplification 
is available in the calculation of the quantum mechanical part of the current J q . 

To begin with, the wavefunction of the quantum mechanical minisuperspace mode (p is not 
gaussian. Moreover, in the tree-level approximation the graph of the probability distribution 



(2J3) is very flat. It does not have good probability peaks and is even unnormalizable. This 
means that the tree-level quantum averages (AQAQ) ree are badly defined. Beyond the 
tree-level approximation the situation can be improved, because they should now be defined 
with the aid of the reduced density matrix 

(AQ(t)AQ(t')) = tr [AQ{t)AQ{t') p] , (3.49) 

p = p(<p, ip') = J df*(<p, /)*V, /), (3-50) 

which originates from tracing the /-variables out and includes loop corrections. As shown 
22 , p4] , 41] , the diagonal element of this density matrix - the distribution function of 



m 



p((p) = p((p,(f), (3.51) 



is given in the approximation of a gaussian integral by the effective action algorithm (|2^ 
Effective action contributes the factor that can generate a sharp probability peak (|2.12 ) with 



the dispersion A defined by (|2.11| ). With this modification the quantum correlators become 
well defined, being expressed in terms of (AipAip) ~ A 2 < oo. Certainly, this improvement is 
achieved by exceeding the precision of the one-loop approximation - badly defined tree-level 
quantum correlators become finite due to one-loop contribution (therefore, in their turn they 
effectively contribute to the radiation currents two- loop quantities). But overstepping the 
conventional rules of the loop expansion is justified here because it reflects the underlying 
physics of the slow roll dynamics. 

Point is that the inflaton field in models satisfying the slow-roll conditions effectively 
represents the massless scalar field - its mass is roughly proportional to the slow-roll smallness 
parameter 0. But massless scalar fields do not have a well defined DeSitter invariant 
vacuum [17]. This fact, in particular, manifests itself in the unnormalizability of the tree- 
level wavefunction exp[=|=/(y?)/2], absence of its local maxima, etc. As we see, loop effects 
render this state a quasi-gaussian nature (|2.12|) and thus justify the improved semiclassical 



expansion. A major part of the paper in what follows deals with the direct calculation of 
the quantum mechanical radiation current J q and its physical implications. 



9 Note that the notion of S e s[Q), as a generator of equations for expectation values, is legitimate only 
within the local derivative expansion. For nonlocal contributions this action does not exist at all - there is no 
mean field functional that could yield by the variational procedure effective equations for expectation values 
|36t . The reason of this is that for nonlocalities the analytic continuation from the Euclidean to Lorentzian 
spacetime is not unique - in addition, it requires setting the retardation boundary conditions for nonlocal 
form factors (see |3(| and cf. the previous footnote). These boundary conditions prohibit the existence of 
the effective action for expectation values. 
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4. Quantum Cauchy problem: tree level approxima- 
tion 



Loop expansion for effective equations is essentially perturbative. Therefore, we solve 
them by iterations starting with the classical solution. Then, in the one-loop approxima- 
tion the radiation current can be calculated on the classical background - the lowest order 
approximation for the mean field. Here we pose the initial conditions for this solution that 
follow from the no-boundary and tunneling cosmological wavefunctions. 

In this and the next section we work with the model of minimally coupled inflaton field 
having a generic potential V(<p) (we reserve the notation as opposed to the notation ip 
for the non-minimal inflaton). This general framework of the Cauchy problem for the cos- 
mological background and perturbations can be easily extended to include the non-minimal 
model by reparametrizing the latter to the Einstein frame fTI , |38| , and this will be done in 
Sect. 6. Thus, we begin with the action 



d 4 x g 



1/2 



(4.1) 



Under the (unperturbed) ansatz for spatially homogeneous metric (|3.22 ), it takes the min- 
isuperspace form 



S[a,cf>,N}: 



dtNa 3 ^y 



K^CL 



a 2 



N 2 a 2 



+ - 



K 



57T 
nip 



Classical equations for this action in the cosmic time gauge, N = 1, read 



(4.2) 
(4.3) 



5S 



7 SN 



3 

K 



- \ - v W 



1 5S 
3Na 2 ^f ~5a 



1 

K 



1 da 
^r + 2- + - 



+ y - V(4>) 



0. 



(4.4) 
(4.5) 

(4.6) 



The first of the eqs. (|5.1| ) represents the nondynamical Hamiltonian constraint. In terms 
of the momenta conjugated to a and 0, ILj, = y / 7a 3 0/iV, II a = —Q^yad/nN, this constraint 
has the following form 



H(a, 



n a ,n^ 



H 



-II 2 + 



12a^ a 2a 3 ^f 



nj + a 3 v/7 



V(<f>)- 



0. 



(4.7) 



which at the quantum level in the coordinate representation of the quantum minisuperspace, 
II a = d/ida, = d/id(p, gives rise to the minisuperspace Wheeler-DeWitt equation on 



H(a, 0, d/ida, d/id(j>)^f{(j>, a) = 0. 



(4.{ 



There are two well known semiclassical solutions of this equation - the so-called no- 
boundary and tunneling wavefunctions. In the approximation of the inflationary slow roll 
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(when the derivatives with respect to 
to a) these two solutions read |2F 



are much smaller than the derivatives with respect 



C™(a 2 ff 2 (,«-l)- 1 / 4 exp 



- 2 m 



COS 



S(a, 



+ 



71 



C T (a 2 H 2 (<P) 



-1/4 



exp 



+-I{<f>)+iS{a, 



+ 



ITT 



(4.9) 
(4.10) 



They describe two types of the nucleation of the Lorentzian quasi-DeSitter spacetime (de- 
scribed by the Hamilton- Jacobi function S(4>, a)) from the gravitational semi-instanton - the 
Euclidean signature hemisphere bearing the Euclidean gravitational action J(0)/2 



S{<f>,a) 



7im 2 P 



Tim, 



2H 2 ((j)) 



[a 2 H\<f>) - 1) 



3/2 



(4.11) 



The size of this hemisphere - its inverse radius - as well as the curvature of the quasi- 
DeSitter spacetime are determined by the effective Hubble constant, a/ a ~ H((f>), driving 
the inflationary dynamics of the model 



W(<p) KV((f>) 



3m 2 P 



(4.12) 



The nucleation of the Lorentzian spacetime from the Euclidean hemisphere takes place 
at a = 1/H(<f)). This domain forms the one- dimensional curve in the two-dimensional su- 
perspace. Its embedding equation can be written in the form 



X(0,a) = a 



kV(</>) 



0. 



(4.13) 



The dimensionality of this subspace coincides with the number of physical degrees of freedom 
in the minisuperspace sector of the model. The intrinsic coordinate on this subspace becomes 
the physical coordinate and the restriction of the Dirac wavefunction \I/(0, a) to this subspace 
becomes the physical wavefunction, provided one takes care of a proper relation between the 
quantum measures on the original superspace and the physical subspace. For a generic 
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constrained system, the details of such a quantum reduction can be found in [^4[ 32, 
Here we just briefly repeat it for our model. 

Let us identify x(0, a) in ( 4.13|) with the gauge condition fixing the time reparametrization 
invariance in the theory (|4.2j) and choose <j) as the physical coordinate. Then, according to the 



formalism of | 32]j , the physical wavefunction *f?(4>) in the one- loop (linear in h approximation) 
can be obtained from the semiclassical Dirac wavefunction ^(0, a) by the transformation 



^(4 > ) = \{ X ,H}\^(4 > ,a: 



x(<M=o 



(4.14) 



Here we distinguish the original Dirac wavefunction in 2-dimensional minisuperspace from 
the physical wavefunction by the number of their arguments. The factor \{x, H}\ - the 
Poisson bracket of the gauge condition with the first class Hamiltonian constraint - is the 
Faddeev-Popov determinant which should be calculated at the semiclassical values of mo- 
menta, n a = daS 1 ^, a), n^, = d^S((p,a). In the slow roll approximation, when the n^- 
momentum is negligible, this factor equals \{x, H}] 1 ^ 2 ~ (H 2 a 2 — l) 1 / 4 and, thus, cancels 
the preexponential factors in eqs. (|4.9|) - (|4. 10|) divergent at the nucleation surface (^4.13|) . 



16 



Thus, the physical wavefunction on the nucleation surface ( }4.13|) which should be regarded 
as the Cauchy surface in minisuperspace reads as 



*jvb,t(0) = CjvB.rexp 



t\h.4>) 



(4.15) 



minus and plus signs related respectively to the no-boundary and tunneling states. It is well 
known that the graphs of these wavefunctions are very flat for the situations when the slow 
roll approximation holds (equivalent to small (^-derivatives). Therefore, they are generally 
not normalizable and do not have good probability peaks that could be interpreted as a 
source of initial conditions for inflation. The inclusion of loop terms via eq. (|2.4j) might lead 
to the normalizability of the wavefunction and, for the model of the nonminimally coupled 
inflaton, even yield a sharp probability peak of the above type. Then, the expectation value 
of the inflaton <fi = ((f)) becomes finite. It is determined by the location of this peak and 
serves as the initial condition for the classical extremal that will be used as the background 
for the calculation of the one-loop radiation currents. 

The second initial condition for this classical extremal - the time derivative of the inflaton 
- arises from the expectation value of the physical momentum conjugated to 0, p^ = (p^>). 
In view of reality of the initial density matrix (|3.50|) this expectation value is vanishing 

-fd*\±**,40l_ = 0. (4.16) 

From the Hamiltonian reduction of the symplectic form in the gauge a) = it follows 
that the physical momentum expresses in terms of the original momenta 



Tl a da + H^dcj) - 

p^ = ^ - UaX^/Xa, X<f> = 9 <l>X, Xa = d a X- (4.17) 

Therefore, for p^ = 0, IT^, homogeneously expresses in terms of Il a and, after plugging this 
relation into the Hamiltonian constraint ( f4.7| ), it implies that at the initial Cauchy surface 
11^ = and Il a = 0. Thus, the full set of initial conditions for the classical background reads 



5. Cauchy problem for cosmo logical perturbations 

In this section we pose the Cauchy problem for quantum cosmological perturbations 
propagating on the classical background of the previous section. First, the set of pertur- 
bations is reduced by the technique of 2(J to the set of linearized invariants of spacetime 
diffeomorphisms, and their quadratic action is constructed. The ghost nature of their min- 
isuperspace sector is revealed and the original perturbations are built in terms of invariants in 
the Newton gauge. Then, quantum initial conditions for perturbations are obtained with the 
aid of the linearized version of the minisuperspace gauge introduced above. Again, we con- 
sider the minimal model which will be later, in Sect. 6, reparametrized to the non-minimal 
curvature coupling. 
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5.1. Hamiltonian reduction to the physical sector 

Here we start with the physical reduction for cosmological perturbations on the classical 
background of Sect. 4. In the main, we follow the notations of [E0| where this reduction 
was presented in much detail. In particular, we use the conformal time denoted by rj corre- 
sponding to iV = a(r]). In this gauge the classical equations of motion ([4.4|)-(fL6"D have the 
form 



-H 



n 2 



2 + K 

i -n' = 



a 2 V = 0, 



J2 



b" + 2H<p' + a% = 0, 



(5.1] 



where primes denote the derivatives with respect to the conformal time, subscript implies 
the partial derivative with respect to the inflaton, = d^V^tfi), and 7i is the "conformal" 
Hubble constant 



n = 



a 



(In 
drf 



(5.2) 



related to the Hubble constant in cosmic time H by the equation 7i = aH. 

The cosmological perturbations (hij, A, Si, 5(p) of metric and inflaton field are introduced 
according to the ansatz 



d<? 

""'total 



-(1 + 2A)drf + 2S i dx i dr] + (7^ + h ij )dx i dx j 



0total = + 50, 

hij = -2ipjij + 2E\ij + 2F { i\ j - ) + tij, 
Si = ViB + Vu V 4 F'' = ViV 1 = ti 



(5.3) 

(5.4) 
(5.5) 



0. 



They consist of the scalar perturbations (if),6<j), E, A, B), transverse vector perturbations 
(Fi, Vi) and transverse-traceless tensor ones Here Vj denotes the spatial covariant deriva- 
tive. 

Spatially homogeneous modes from the minisuperspace sector, upon which we focuse in 
this paper, belong to scalar perturbations. As discussed above, the inhomogeneous modes 
which contribute to the J* radiation current can be treated by the effective action method 
and do not require a manifest physical reduction. Thus we consider only the scalar sector. 
After constructing the quadratic part of the action in terms of scalar perturbations one 
introduces the momenta conjugated to (ip, E) 



n<50 

He 



2a 2 ^7 



-3{i>' --(/>' 6(f> + HA 



A{B - E') 



a 2 ^(64>'-(f>'A), 



2a^/j 



A 



+ HA-(B- E') 



(5.6) 



and finds out that (A, B) play the role of Lagrange multipliers to the linearized Hamiltonian 
and momentum constraints 



C A = -HTlip + + a 2 ^ 

Cb = n^, 



Dip + W - <f)") 

K 



(5.7) 
(5.8) 
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where D is the following modified covariant Laplacian acting on a closed 3-sphere with the 
metric 7^ 



D = A + 3, A = Y j V i V j 



(5.9) 



The constraints ( p.7| ) generate the diffeomorphisms in the scalar perturbation sector with 
respect to the vector-field parameter A 1 ' = (A , V*A) 



5 x ip = 



-n\°, S X (S4>) = </>'A°, 5 X E = A, 

^H D \^ s x U S4l = a 2 ^((j)" - U(j)')\ , 5 x Ii E = Q 



(5.10) 
(5.11) 



accompanied by the transformations of the Lagrange multipliers 5 X A = (A )' + HX°, 5 X B - 
A' — A . There are two obvious invariants of the gauge canonical transformations (|5.10|) - (|5.11 



n, 



K(j)' ^ 



(5.12) 
(5.13) 



It turns out that after solving the constraints, Ca = Cb = 0, ( |5.7|) - (|5.8|) , with respect to LT^ 
and He and feeding the result into the canonical action the latter entirely expresses in terms 
of these two invariants. Moreover, they play the role of a single pair of canonically conjugated 
variables in the physical sector pPf : on the constraint surface in phase space the original 
symplectic form goes over into the physical one, IL^ip' + Hs<t>8(j>' + IIe-E" = 11^ St' + (...)'. The 
corresponding canonical action quadratic in (\l>, H^) reads 



drj 



2a 2 V7 



.-2^/2 



2a 2 Ji % 



aV7 

K 



K 1 

—t, — — n^— 

4a 2 v 7 ^ D 



(5.14) 



Somewhat simpler form this action acquires in terms of the new variables {q,p) related to 
(|5.12|) - (|5.13 ) by the canonical transformation 



la T HI 
* H 

K(p cp' dyJ^D 



p - — />»w —11, (5.15) 
In terms of them the quadratic action in the physical sector of scalar perturbations looks as 



S[q,p] 



drj 



(j)" K<f) 



/2 s 



pq +pq \^ + 

K(j)' 2 < 



AH 



1 1 
+ —^p—p 



8H 2 



1 - 



K(j)' 



/2 X 



4 



ViqDq - -y/j(Dqf 



(5.16) 



With the extremal expression for the momentum 



P 



70' + K<j) u /m)q 



(5.17) 



19 



the Lagrangian form of this action is even shorter 



S[q]\ 2 = lfdnyfy(-Dq) 



d 2 „/l\" Kef)' 2 



drf T V y 2 



(5.1* 



The invariant field q here is well known from the theory of cosmological perturbations pOf . 
It is actually given by the so called Bardeen variable |39], |2"0|| , &h — 7~t{B — E') — if), q — 
-2a$> H /K(j)'. 

Note that the operator D given by ( |5.9| ) is negative definite except for two modes: the 
zero mode corresponding to the Laplacian eigenvalue A = —3 and the spatially homogeneous 
mode for A = 0, D = +3. In view of the overall factor — D the zero mode does not enter the 
action at all, while the homogeneous mode enters ( |5.18| ) with a wrong sign - its kinetic term 
is negative. Thus, this is a ghost variable signifying the classical instability of the model. 
This instability at the linearized level is nothing but the manifestation of the inflation which 
is a huge instability phenomenon incorporating the runaway modes. In contrast with a 
conventional wisdom of the S-matrix theory, this instability should not be regarded as an 
irrecoverable flaw of the theory, because we know a nonlinear damping mechanism that 
provides an exit from the inflation stage in case of the inflaton field rolling down to smaller 
values of the potential. In particular, no special measures like introducing the indefinite 
metric should be undertaken to eradicate this phenomenon. Homogeneous fluctuations of 
the inflaton field do not have a particle nature and one should not take care of guaranteeing 
the energy positivity of their excitations. Therefore, this mode can and should be quantized 
in the coordinate representation with positive metric in the Hilbert space. 

A single spatially homogeneous mode q{r]) contained in the full set of 

q(x) = (q(r)), q(r),x)) (5.19) 



corresponds to the D = +3 eigenvalue of the operator ( |5.9|) in the action ( |5.18|) . It also 



satisfies all the above relations with a simple ultralocal substitution D = +3 and is actually 
responsible for the perturbations in the minisuperspace sector of the cosmological model. 
Indeed, from the metric ansatz (|5.3|) , (|5.5|) it follows that spatially homogeneous variables 
ip(r]), 54>(r]) and A{rj) induced by qirf) generate the variations of the scale factor, inflaton 
field and lapse function 

5a= -aiP + 0(ip 2 ), 5(f), 5N = aA + 0(A 2 ). (5.20) 

Actual expression for ip(f]), o~4>{t]) and A(r]) in terms of {q{f])iP{ff}) depend on the par- 
ticular gauge chosen for minisuperspace variables. In what follows we will need two types of 
such gauges. One will be used for gauge fixing the dynamical evolution of perturbations as a 
function of dynamically evolving invariant variable q(i]). Another gauge serves as a part of 
the quantum Cauchy problem - as shown in the previous section, it facilitates the quantum 
reduction to the physical sector and relates the wavefunction to the initial conditions for 
both the classical background and the homogeneous perturbation variable q(i]). The first 
gauge may coincide with the second one. However, its use is strongly biased by practical 



necessities of the theory of cosmological perturbations {40[ and, therefore, is usually chosen 
to be the Newton gauge which is essentially different from the minisuperspace gauge of Sect. 
3. Thus we consider these two gauges separately. 



5.2. Newton gauge 

Newton gauge is widely used in the theory of cosmological perturbations to express 
them in terms of the Bardeen invariant q. The Newton gauge for spatially inhomogeneous 
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modes reads 



B = 0, E = 0. {5.21) 



From the equations for momenta ( |5.6|) and the momentum constraint Cb = this implies 
that 

P. 2 I — 

U^ = -?!ivlty-l < f ) 'S < f> + HA) = 0. (5.22) 

In the spatially homogeneous sector of the theory, where the contribution of B and E is 
missing (they enter only differentiated with respect to spatial coordinates) the latter equation 
should be regarded as the definition of the Newton gauge. This gauge involves only the phase 
space coordinate - the momentum IL^ - and, therefore, it is unitary and can be identified 
with the gauge ( |3.13|) . 

Canonical equations of motion for (^,11^), which follow from the action (|5.14j) , in this 



gauge have a simple corollary (atp)' /a — K(f) > 5(f)/2 = 0. When compared with fl5.6|) , 11^ = 0, 
this corollary yields the main relation in the Newton gauge 

A — ip. (5.23) 

Then one easily expresses all the perturbations in terms of the physical phase space variables 
(q,p). On substituting the Lagrangian value of the momentum ( |5.17|) these expressions 
finally simplify to 

K(b' 

*l> = ir-Q, (5-24) 
"-<*«>' (5.25) 



a<p' 



Eqs. ( |5.23|) -( |5~25"1) form the needed set of relations ( p.!6|) of the physical reduction for 
minisuperspace perturbations. 

5.3. Minisuperspace gauge 

The minisuperspace gauge x(o> 0) = of Sect. 3 was used for the physical reduction of 
the minisuperspace wavefunction and for establishing the tree-level initial conditions - for 
the classical background. Let us now use it in order to find the initial conditions for (q,p). 



The linearized minisuperspace gauge condition (the gauge ( |3.13| ) in condensed notations of 
Sect. 3.1) gives the perturbation of a in terms of 5<f). Taking into account the relations 
( |5.20 ), expressing Sa in terms of the perturbation ip, we get 



$ = (5.26) 

The corresponding reduction of the symplectic form gives the expression for the physical 
momentum ps^ conjugated to 5 (ft 

U^' + U s< p8(j)' = p5<j,S(j)' + p s<t> = H S( p + 11^^- . (5.27) 
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Then, by solving the linearized constraint Ca = one easily finds 11^ and 11^ as functions 
of (5<f), ps(/>) and, via the formalism above, proceeds to the final transformation relating (q,p) 
to 



q{v) 



i 

M 



(H 

\ K(f)' 



"\a 2 Xa H 



s<Kv) + 



(5.28) 



v 4M 
M = (p'x^ + Haxa 



+ 



2X 



2 1 



(5.30) 



Note that this relation is written down in the homogeneous sector which is emphasized by 
the time arguments of the phase space variables. In the right hand side of these equations 
the spatial homogeneity manifests itself in the particular value of the operator D, D — 3. 
One can easily check that this transformation is canonical, {q,p} = {5(f>,ps^} = 1, and 
invertible. Inverting it, one can find all the minisuperspace perturbations as functions of 
q and q' , similarly to the relations ( |5.23| )-( |5.25| ) in the Newton gauge. However, the goal 
of working in the gauge ( 4.13 ) is somewhat different. We shall need the relations ( [5.28 )- 
( |5.30| ) in order to express initial conditions for (q(rj),p(rfj) in terms of initial conditions for 
(5(f>, psj,)- The latter in their turn follow from the cosmological wavefunction in the physical 
sector. 

We begin by noting that at the initial moment of time the following relations hold 



n = r ] + 0(r ] 2 ), (f)' 



-r] + 0(r] ), — 



kV 



Mi 



kV 
3l4 : 



/ 3 \l/2 



0. 



Using these relations in eqs( |5.28j )-( |5.29| ) we obtain the asymptotic behavior of the invariant 
variables (q, p) for r] — > in terms of the physical variables of the minisuperspace gauge 
fixing (5(f)(rj),p 5 <p(ri)) 



IT 



1/2 



1/2 



P 



1 rj. 

— 5d> H 

3r] Y 9 1 

9 



9e 2 /16 -/y' 



kV 



;i-9e 2 /32)^, v -^0. 



(5.31) 
(5.32) 



Just to emphasize the role of the slow roll expansion we retained here the corrections pro- 
portional to the smallness parameter 



8 VI 



3kV< 



<C 1. 



(5.33) 



In what follows we shall systematically discard such corrections retaining only the leading 
order of the slow roll expansion. 

Important peculiarity of the behavior (|5.31|) -( |5.32|) is its singularity at rj — >■ 0. This sin- 
gularity is, however, an artifact of the definition of the invariant variables ( |5.15|) nonanalytic 
at </>'—>• 0, rather than the manifestation of some physical inconsistencies. To see it, one 
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can decompose the general classical solution for q(rj) in the sum of two linearly independent 
solutions of the equation of motion for the action (|5.18 ) 

q(v) = c+q + (v) + c-q_(ri), (5.34) 

Because of <f>'(l/(j)')" ~ 2/r/ 2 , r\ — > 0, the initial moment rj = is a singular point of this 
differential equation, at which one of the two solutions, q_(rj), diverges as 1/0'. One can 
make a singular rescaling, 

g=5 (5.36) 
to a new variable Q(?y) which is finite at this point. It satisfies the equation 

+ 2^-— + ^- + 3 Q(r/) = 0, (5.37) 



y dr] 2 (J)' dr] 2 
and has as two solutions the following regular functions 

Q + {r,) = rf(l + 0{rf)), (5.38) 
Q_(rj) = l-3r] 2 /2 + 0(r] 3 ). (5.39) 

Substituting the decomposition ( |5.34| ), with q± related to Q ± by ( |5.36| ), to the left hand 
sides of ( [5.31| )-( |5"l?2j ) one obtains the system of equations for c± with singular coefficients. 
This system, however, has a regular solution in terms of the initial conditions for physical 
variables (50(0), p^(0)) 

This basic relation will be used throughout the rest of the paper to express the Heisenberg 
operators of quantum perturbations AQphys(^) and AQ(rj) in terms of the Schrodinger op- 
erators, 5(f)(0) = 5<fi, p6cf>(0) = d/id(S<p), and then find the quantum avera ges of their bilinear 
combinations in the inflaton representation of the initial density matrix (|3.50|) . 



6. Non- minimal model 

In what follows we go over to the model ( I 2.5 ) that has a good peak-like behavior of 
the initial distribution function of the inflaton 
action in this model can be rewritten in the form 



[Tql . The inflaton-graviton sector of the 



S[g^ <p] = J d'xg 1 ' 2 {-Vfo,) + U{fp)R - ^(^(V^) 2 } . (6.1) 

In fact, the curvature (and derivative) expansion of any low-energy effective graviton-scalar 
action can be truncated to this form with some coefficient functions of the zeroth and first 
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order in the curvature - the scalar field potential V(<p), the effective (^-dependent Planck 
"mass" !6irU(ip) and the coefficient of the inflaton kinetic term G(<p). In the classical model 
Q2.5|) these functions have a particular form 

U^) = ^- + \\i\v\ (6-2) 

l07T 2 

v m = — y~ + -4-: (6-3) 

G{tp) = 1. (6.4) 

It is well known that the action ( |6.1| ) can be transformed to the Einstein frame by a special 
conformal transformation and reparametrization of the inflaton field (g^jCp) — > (g^,^), 

S[g, v M = J d^xt 12 + ^R(g»») - ^(Vy?) 2 } • (6.5) 

In what follows, we shall denote the fields and other objects in the Einstein frame of the 
non-minimal model by bars and identify them with those of the minimal model considered 
in Sects. 3-5. In this way we reduce all the calculations, Cauchy data setting, gauge fixing, 
reduction to the physical sector, etc. to the algorithms derived above for the case of the 
minimal model. 



6.1. Reparametrization to the minimal frame 

The transformations relating the actions S[g^ u , ip] and S[g^ u , (p\ are implicitly given by 
equations [B8L ITTJ 



l6nU(ip) 



m. 



'dtp\ 2 
dip) 



m 2 P GU + 3U 2 



16tt 



U 2 



(6.6) 
(6.7) 



where, similarly to previous sections, ^-subscripts denote the derivatives of the coefficient 
functions with respect to the inflaton, V v = dV/d<p, V w = d 2 V/dip 2 , etc. The action in 
terms of new fields ( |6.5| ) has a minimal coupling and the new inflaton potential 



V{0) 



m, 



16tt) U 2 (ip) 



(6- 



<p=ip(<p) 



For the coefficient functions (|6.2| )- (|6T3| ) the explicit reparametrization between the frames 
can be found for large value of the nonminimal coupling constant |£| 3> 1 and small value 
of the parameter m 2 P / |£|y? 2 <C 1 |JTTJ 



V{0) 



mp 

W 2 



exp 



Xmp 



256tt 2 |£| 2 



1 + 5 m 2 



2 
p 



An \£\ip 2 



+ ... 



v=v(v) 



(6.9) 
(6.10) 
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where we have retained only the first order term in mp/|£| f 2 . In view of ( |6.9|) , for large if 
this potential exponentially approaches a constant and satisfies the slow roll approximation 
with the expansion parameter [pH 



-1/2 



3tt ^) 



1 + 



ml 



eiei. 



(6.11) 



which justifies the smallness of the parameter mp/|£|y? 2 chosen above. 

Let us now consider the minimal model of Sects. 3 - 5 as an Einstein frame of the non- 
minimal model and label all the objects of the minimal model - the metric, inflaton field, 
scale factor, conformal time, cosmological perturbations and the minisuperspace gauge fixing 
conditions - by bars 



9w, <P = Pi a, rj, 

8<p = 8f, A, n^, fity, x(a,<p), 



(6.12) 



as opposed to the objects in the original - non-minimal - frame: g^ u , if, a, rj, ip, 8f, A, x(a, if). 
Comparing the metrics in these frames, perturbed by the cosmological disturbances from the 
scalar sector, 



ds 2 
ds 2 

ds 2 



[1 + 2A)drf 



2ip)^ijdx l dx 3 



-(1 + 2A)dff + (1 - 2ip)^ ij dx i dx j 



16ttU 



ds 2 



mi 



one finds the relations between these two sets of variables 



\ 



16ttE%) 



ml 



\ 



ml 



T] = T], 



= 8 if 



m 2 P U + 3U' 2 

167T U 2 




If = If 



A 



U, 



8f 

ip , 

<P 

A + 5 ^ 



(6.13) 
(6.14) 

(6.15) 

(6.16) 
(6.17) 
(6.18) 



where the last three relations hold in the linear order of perturbation theory in cosmolog- 
ical disturbances. The canonical momenta Us v and LT^ obviously transform by the rule 
contragradient to (|6.16|) and ( |6.17|) . 

We also need the frame transformation between the physical sectors defined in the min- 
isuperspace gauge of Sect. 5.3. The gauge condition itself transforms as a scalar - only in 
this case it represents one and the same Cauchy surface, written in two different coordinate 
systems on minisuperspace 



(6.19) 



As regards the reparametrization beteween these coordinate systems, ( |6.9| ) and ( 6.15|) , it has 
a general form 



f = f(if), 



a{f,a) 



(6.20) 
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mixing the inflaton and the scale factor only in the transformation of a. Therefore, the 
linearized perturbations of a and (p and their momenta in both frames are related by a 
triangular transformation 



5if = ^p 5ip ' ^ 



1 da a da 

— ^~ 5l P + - o- ^> 
a Off a da 



n 



Sip 



1 da 



— — — n^, 

dip a dip 



a da 
a da 



n 



(6.21) 
(6.22) 



The physical momentum ps v expresses in terms of phase space momenta by the barred 
version of eq. (|5.27 ). Then, in view of the above relations, one easily finds 



dip 



(6.23) 



This equation holds exactly for an arbitrary choice of the gauge condition function x( a , v), 
and this is a corollary of the triangular form of the transformation ( |6.20| ). In our nonminimal 
model with |£| 3> 1 this implies the following simple relation between the physical sectors in 
two frames 



Sip 




3 dip 
m P — , p Stp 




4tt ip 



3 mp 



(6.24) 



6.2. Quadratic order currents 

In the minisuperspace sector of the non-minimal model Q 
tional derivatives of the classical action read 



(N(t),a(t),<p(t)) the func- 



' II = 6t%) ( i + + GUM*- " £ " V{<P), 



a 3 ^7 SN 
1 SS 



Na 3 -^ Sip 



6-ip 

a 



a 



a 



a 2 " 



* i I I a 
6U v (<p) - + - + - 



a 2 a a 2 



Vp{<p), 



(6.25) 
(6.26) 



where dots are used to denote the parametrization invariant derivative a = da/Ndt, ip 
dip/Ndt. Now we use the perturbed ansatz ( |6.13| ) for total minisuperspace variables in these 



equations, N 



N 2 

iv tot 



a 2 (l + 2 A), a 2 



not 



a 2 (l - 2-0), ip — >■ ip tot = ip + Sip, and 



carefully expand the first order variations of the classical action up to the second order in 
perturbations (A, ip, Sip) on the classical background. The result reads as follows 



1 
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In 



A 2 (24UH 2 + 2AU tp Hip' - 2ip 
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V) 
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+ipSip 
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6S_ 

Sip 



+6U tpip H5cpScp' - 



-- A 2 (2a 2 V 9 - 12U V ) + AA' (18U V H - &p') 
+A^j (6a 2 V v - 24^) + (Aip)' (18U V H - 3<p') + 6£^(A0')' 



+A5tp 



6t/ w (l 



H 2 - H') - aV w + 2HA8ip' 



-QU w HA'5<p + (A5tp')' 

-\2u^ 2 + w (i&u v n - &p') + eUyW 1 
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3a 2 V 



mj w ni)'bip 



3U. 



(l + H 2 + H') - -aV v 



(6.27) 



(6.28) 



Here the unlabelled variables correspond to the minisuperspace background in the conformal 
time gauge, N — a, t — rj, primes denote derivatives with respect to rj and TC = a' /a 
denotes the conformal time background Hubble constant. The background satisfies classical 
equations of motion which were used to simplify the coefficients of the above quadratic forms 
in AQ = (A,ip,5(p). 



6.3. Quantum rolling force: effective action and minisuperspace 
contributions 



In the presence of the spatial densities of one-loop radiation currents 



JN 



a 3 Jl 



In 

5S 

5ip 



(6.29) 
(6.30) 



N and ip components of the effective equations of motion in the non-minimal model read 

;2\ • -2 



. a ip 
a ~2~ 



a a 



-<p - 3-<p + m v I — + - + — I - V v + j v = 0. 
a \a z a a z ' 



(6.31) 
(6.32) 



In view of eqs. (|3.10| ) and ( |3.20| ) their a-component expresses in terms of the above two 
ones, so that eqs.( p.31 )-( |6.32j ) in a consistent manner exhaust the quantum dynamics of the 
mean fields. Differentiating the first of them with respect to time one obtains the system 
of two equations for a and (p. Substituting the solution of this system for a into the second 
equation one finally has the equation for the mean inflaton field with quantum contributions 
to the friction term and the rolling force 



(6.33) 
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H<P, a, 0) = 2V Hl 3^ - \p 2 ^ HU + K) + F loop (<p, «, 0, A), (6.34) 
F loop (^, a, 0, a) = Z7q ^ ^ - 2U,j N - • (6.35) 



The first two terms in eq.( 6.34 ) represent the classical rolling force, the (p 1 contribution 
belonging to the subleading order of the slow roll expansion. As regards the quantum part, its 
radiation currents in the one-loop approximation split into the contributions of the quantum 
mechanical sector and the field sector of spatially inhomogeneous modes, ji-ipop = 3 q + 3 



(cf. eq.( |3.34|) ). According to the discussion of Sect. 3.1, see eqs. fl3.45|) - (|3.48|) , the /-part 



of the current can be absorbed by the replacement of the original classical action with 
the effective one ( |3.48|) . This implies the replacement of the classical coefficient functions 
V((p),U((p),G(ip), ( |6.2|) - (|6.4|) , by their effective counterparts 

S^M = J d'xg 1 ' 2 {-V cS &) + W*(<p)R- l -G^{ V ){V V f + ...} , (6.36) 

and truncation of the (generally infinite) series to the first three terms. This truncation is 
based on two assumptions - the smallness of inflaton derivatives due to the slow roll regime 
and smallness of R/m 2 ^ rt - the curvature to particle mass squared ratioQ. Thus, with this 
approximation, the effective equations of motion in our non-minimal model take the form of 
(|6.31 ) and ( |6.33|) with V e e((p), U e s(ip), G e ff(y?) replacing V(<p), U(<p), G((p) and the radiation 



currents Jtv, j v saturated by the contribution of the quantum mechanical mode, jjv, j^. The 
resulting rolling force in the leading order of the slow roll expansion becomes the sum of the 
force induced by the effective action, F eS , and the quantum mechanical force, F q , 



F = F + F\ (6.37) 

9 T/eff j reff 7"/"effT/cff 




^ = W^m-^ r A-^t). (6-39) 

The one-loop calculation of V e s((p),U e s(ip) and G eS (<p) and the effect of F eS on the 
inflationary dynamics have been studied in [|IJ. This effect is qualitatively different for the 
no-boundary and tunneling cases and briefly looks as follows. For the no-boundary state the 
one-loop corrections in the distribution function add up to form the full Euclidean effective 
action 

Pnb (v 9 ) = const ex P [- r (f)} , (6 AO) 

VeS{<P) 

Its value on the DeSitter instanton follows from that of the classical Euclidean action, 1(f) = 
— 96ti 2 U 2 /V , by replacing the classical coefficient functions V(f), U( f) an d G(<p) by the 
effective ones, coinciding with those of the Lorentzian effective action Q6.36 ). Therefore, by 



10 The Schwinger-DeWitt expansion involves the inverse powers of masses of particles of constituent quan- 
tum fields. The latter acquire their masses via the Higgs effect due to the interaction with the inflaton, so 
that this ratio becomes order of magnitude A/|£| <C f Q. 
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the direct inspection of ( |6.38|) one observes that the effective rolling force in the no-boundary 



case is proportional to the derivative of the distribution functionf 1 ] 

1 (V cS ) 2 d 




^1(1 + 6) J, S) +0{m °), (6.42) 



48vr£ 2 r \ <pj 



and, thus, vanishes at the probability peak ipi. The no-boundary peak is realized for 1+5 < 0, 
therefore the point ifj turns out to be an attractor - quantum terms in effective rolling force 
lock the inflaton at its constant initial value and give rise to infinitely long inflationary 
scenario with exactly DeSitter spacetime. 

In the tunneling case, the distribution function is not related to the overall effective 
action, because its tree-level part has a wrong sign. The probability peak exists in the 
opposite range of the parameter (|2.9|), 5 > — 1, and the rolling force 



ff= _X^ + t) J l + S ]+0{l/m (6 . 43) 



has the quantum term which initially doubles the negative classical part. Therefore, the 
inflaton starts slowly decreasing under the influence of this force, and the tunneling state 
generates a finite inflation stage with the estimated e-folding number (|2.15|) . These con- 



clusions disregard the contribution of the quantum mechanical radiation currents, and we 
proceed to their calculation. 



6.4. 



Quantum state of the minisuperspace perturbations and their 
correlators 



The calculation of quantum averages in the quadratic currents (|6.27|) and (|6.28| ) requires 
the set of quantum correlators (|3.49|) of bilinear combinations of minisuperspace disturbances 
AQ = (A, ip, 5ip) and their derivatives. For this purpose we, first, need the reduced density 
matrix of the inflaton field p(ip, if') in the non-minimal model and, second, the expressions for 
the Heisenberg operators AQ(r]) in terms of the Schrodinger operators of initial perturbations 
and their momenta, 5ip = 8<p, ps v = d/id(5f). 

As we know, the diagonal element of the density matrix has a quasi-gaussian behavior 
( |2.12| ), which is, however, insufficient for averaging the operators involving momenta. The 
necessary off-diagonal elements with one-loop contributions of various massive and massless 
fields have been calculated in HIL fi2]| . It was shown that in the model with a big |£| the 
initial density matrix describes practically pure quantum state and expresses in terms of the 
distribution function 



p(<?, v') | t=0 - vW) vWy, if i » i- 



(6.44) 



The explanation of this property jy] is based on the fact that the decoherence factor D(<p, ip') 
by which the initial density matrix differs from Q6.44|) is a function of the arguments m/H(ip) 

11 The explanation of this observation is simple. In the minimal frame the rolling force is given by the 
gradient of the potential, while the logarithm of the distribution function is inverse proportional to it, 
the com bina tion V eS (</?)/[ U eS (y) ] 2 representing the minimal frame potential in terms of the non- minimal 
objects 
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and m/H(ip') for a quantum field of a mass m. For large |£| the masses of particles generated 
by the Higgs effect give rise to big and predominantly ^-independent ratio m/H(ip) ~ \f\£\, 
so that D(ip, ip') ~ 1. For massless fields a similar conclusion can be drawn because for them 
the role of mas s is p layed by the Hubble constant H(cp) of the quasi-DeSitter background. 
In view of (|2.12| ), the effectively pure quantum state, 



NB,T 



(<p) - \[pnbM}p), 



(6.45) 



in the vicinity of the probability maximum, which is located at ipj, can, thus, be approxi- 
mated by the gaussian packet of small quantum width A 



NB/T 



(Sip) = ^ nb ,t(^i + $<p) 



(2tt) 1 /4^A 



cxp 



_ 4A 2 



(6.46) 



The operators of quantum disturbances in the (^-representation, acting on the wavefunc- 
tion of the above type can be found by c ol lecti ng together several sets of equations derived 
above. First, we use the equations ( |6.16|) -( |67L8| ), relating AQ to the Einstein frame pertur- 
bations AQ. Then, we apply the barred version of equations ( |5.23|) -( |5~25"1) to express the 
minisuperspace perturbations AQ in the Newton gauge as functions of the invariants q and 
q' in the minimal frame. Finally, we use the set of equations ( |5.34|) , ( |5.36|) and ( |5.40D -( |5.41|) 
with barred (minimal frame) potential and physical variables to express these invariants in 
terms of 5<p and pg^. The final result looks as follows 
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(6.47) 
(6.48) 
(6.49) 



where the operators c± with the aid of ( |6.24j) read as 




(6.51) 



Now we are ready to find the quantum correlators necessary for the radiation currents. 
We choose a symmetrized combination of disturbances and their conformal time derivatives, 
AQip = {AQ, AQ', AQ"), in the definition of the correlator 



\AQ l AQ 2 ) = 1 / d(8ip)^*(S^)(AQ 1 AQ 2 + AQ 2 AQ 1 )^(8ip), 



(6.52) 



because in the Hermitian operators of quadratic currents (|3.9|) the products of operator 
valued disturbances automatically enter in symmetrized form (in view of the symmetry of 
the 3- vertex function). The further calculation of the correlators and radiation currents is 
straightforward. However, the general answer that involves the basis functions Q±(r)) for 
arbitrary 77 is still very complicated. Therefore, we separately consider the beginning of the 
inflation epoch, rj = 0, and the late stationary stage of inflation. 
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7. The onset of inflation 



At the onset of inflation the basis functions Q±(rj) have a behaviour ( |5.38| )-( |5.39| ). Using 
it in eqs.( |6.47|) - (|6.49| ) with the operators c± defined by (|6.50|) - (|63lD one easily obtains the 
initial equal-time correlators with respect to the gaussian state fl6.46|) . 

In the leading order of the slow roll expansion those correlators that do not involve 
derivatives (potential type correlators) read 



2 . , . . , o , A 1 



288vr 2 |£| 2 f 
Xipi 1 



288vr 2 |£| 2 /' 

while the correlators of conformal time "velocities" (the kinetic type correlators) equal 

(W)„ = -<W)o = 5^/, 

<^>° = 28^ < 7 ' 2 > 

As we see, these two groups of correlators differ by the power of a special parameter / which 
is inverse proportional to the square of quantum dispersion of the inflaton field 

f = _* L.MVi (73 ) 

J 48tt 2 |£| 2 «;A 2 \1Qtt 2 ) \l + S\' K ' J 

Such a dependence on / reflects an obvious fact that the kinetic type correlators (or corre- 
lators of momenta) in the gaussian state of the form ( |6.46| ) are inverse proportional to A 2 , 
and thus grow with A — > 0, while the potential type correlators (correlators of coordinates) 
are proportional to A 2 , and thus decrease with the narrowing of the gaussian peak. 
The calculation of mixed correlators with one or two derivatives of the form 

( AQAQ' ) = ( AQAQ > O(e 2 ), ( AQAQ" ) = ( AQAQ > O(e) (7.4) 

shows that they belong to the subleading order in the slow roll parameter ( 6.11| ). Finally, 
the additional correlators with three derivatives, which arise in the calculation of dj q N /dt in 
the quantum rolling force (|6.39|) , express as 



(<V)oHW)o = - 2 WV> 
(VV>o = -27W 2 > 0) 



o- 



W) = -2W 2 > , Ty^O. (7.5) 

Although they tend to zero in view of TC(t]) — > at rj — ► 0, their contribution to the quantum 
rolling force is nontrivial because in ( |6.35D they are devided by a/ a = Ti/a. 

Let us now go over to the calculation of radiation currents at rj = 0. Within the slow 
roll approximation, mp/|£|y? 2 <C 1, |£| 3> 1, and in view of a particular form of classical 
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coefficient functions V(ip),U((p), the quadratic currents (|6.27|) -( |6T28D are dominated by the 
following expressions involving both the potential and kinetic terms 



j%(0) 



v 

x^ 



-2(^) + -(#^)-^(V) + (^ 2 ) 

4^> + -(^> + J(iV> 
<p z 



-<W> 



(7.6) 
(7.7) 



On using the tables of correlators above, the radiation currents, contributing to the quantum 
rolling force, take the following final form 



3n(0) 

^(0) = Xtf 



Xipj X 



4 96vr 2 |£| 2 
A / 



96tT 2 |£|< 



1 ~ S f J 



a dj N 
2a dt 



(0) = 0. 



(7.8) 
(7.9) 
(7.10) 



These quantities are strongly suppressed as compared to their classical values, —j^ = 
V((fii) ~ Xipj/i, —j v = V^(ipi) ~ Xipj by a very small 

factor A/|£| 2 ~ AT 2 /T 2 ~ 10~ 10 related to the CMBR anisotropy. Their sign crucially 
depends on the magnitude of the parameter /, (|7.3|) , which in our model is likely to be very 
big, / 3> 1. This follows from the estimate N > 60 on the e-folding number (|2.15| ) and the 
value of |£| ~ 10 4 [p6fl . In this case, the terms proportional to / ~ 1/A 2 , generated by the 



kinetic terms of the radiation currents, (AQ'AQ'), dominate and, in particular, result in 



e q (0) = -j N (0) 



A 2 |l + £| 

isiei 3 
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167T 
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< mp. 



(7.11) 



Interestingly, the sign of the quantum rolling force due to the homogeneous mode is inde- 
pendent of the magnitude of /, because in the leading order of the slow- roll expansion the 

contributions of potential terms, (AQAQ) ~ 1//, completely cancel out 



F q (0) 



U 



36 96tt 2 |£| 



r/>0. 



(7.12) 



In view of the expressions for / and ipi, ( [2.10|) , the magnitude of this force is again much 
smaller than its classical counterpart 
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(7.13) 



Therefore, for the tunneling state it gives a negligible contribution to the effective force 
( |6.43|) . For the no-boundary state, the initial effective force (|6.42j ) vanishes, but the only 
effect that the positive F q (0) can produce in this case is that it shifts the equilibrium point 
from if j to slightly higher value of the inflaton F q (0) + F^ B (ip^) = 0, at which again the 
system will undergo endless inflation. 
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8. Late stage of inflation 



At late stationary stage of inflation the dynamics of the classical background can be 
approximated by the ansatz 

cosh[H (<p)t\, ip ~ tpi, (8.1) 



Hitp) 

with the Hubble constant H((p) approximately linear in (p. In the Einstein frame, it looks 
similar with the Hubble constant which is practically independent of the inflaton H 2 (ip) = 
8ttV (<p) 1 3m 2 p ~ Amp/967r|^| 2 , the cosmic time parameters being related in both frames by 



t ~ t^87r\^\(f/mp. 

The transition period between the onset of inflation and its steady stage can be described 
by solving the inflaton equation with the approximately constant rolling force F and the 
friction term based on the ansatz (|8.1|) for a/ a 

(p + 3if tanh(Ft)^ - F = 0, 

V (0) = Vl , 0(0) = 0. (8.3) 

For late times, Ht 3> 1 (but not so late that the inflaton field evolves too far from its initial 
value), the exact solution to this equation 

V(t) = Vi + ^f 2 ln(cosh Ht) + tanh 2 (#t) (8.4) 
reads as an almost linear function of t 

¥>(«) = <Pi + £gt+ 3^(1 - ln2) + O (e- 2 ^) . (8.5) 

This behaviour corresponds to neglecting the (p term in the inflaton equation of motion 
and solving it for (p, (p ~ F/3H. In our model with the classical rolling force ( |2.14j ), 



ip ~ —AipHe/3 Hip with e ~ m 2 p/\!i\<p 2 ~ the slow roll smallness parameter (|6.11 ). Thus, 



in the lowest order of the slow roll approximation the inflaton field remains constant. 

Let us study the behaviour of the basis functions Q±(t)) for Ht ^> 1. To begin with, 
note that for late times corresponding to exponentially large scale factor the potential terms 
in the wave equation for Q±(r}), ( |5.37| ), can be discarded. The first one, K<p' 2 /2 = 0(e 2 ), is 
small in view of the slow roll regime and the second one, the spatial curvature term 3, is 
small compared to the kinetic terms growing with a, d 2 /dr] 2 ~ a 2 d 2 /dt 2 . Therefore, at late 
times this equation simplifies to 

and has two explicit solutions Q± 

Q-('l) = 1, (8.7) 
QM = N + £ *> N + = ( W 2 (8.8) 
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compatible with the expansions ( |5.38|)-(d.39|) at early times 77 — > 0. (The mismatch between 
the constant function (|8.7|) and eg. ( |5.39 ) has a simple explanation: the term — 3r/ 2 /2 in Q 
of eq. (|5.39|) is induced by the curvature term which we discard at late times, while the 0(rf) 
corrections are due to the perturbation K,ip' 2 /2.) In terms of the cosmic time in the original 
frame, Q + represents a well known growing mode jE| which for late times in the slow roll 
approximation reads as 

Q + = -smhHt. (8.9) 

Thus at late times the operator of the invariant cosmological perturbation is dominated by 
the growing mode 



Q(t) = c+Q + (t) + c_Q_(t) ~ c + Q + (t), Ht > 1, (8.10) 

and in the Newton gauge all the minisuperspace perturbations express in terms of one 
operator c + defined by the momentum ps v , Q6.50Q , 

A = ^ kQ 



2a(t) 




Since 5(p contains ip in the denominator, all the other perturbations in the minimal frame 
are much smaller in magnitude, (A, x/j) ~ 0(e)5ip/mp -C 5(p/mp. Therefore, in view of 
eqs. ( |6. 16|) -( |6. 18 ) the perturbations in the non-minimal frame read 




A~-ip=~\ 5<p~\ —Sip, Hf>l. (8.12) 

V 3 mp V 3 mp 

Another important property of the perturbations in both frames is that to the leading 
order in slow roll they are constant in time for Ht ^> 1. This follows from eqs. (|8.11|) contain- 
ing the exponentially growing functions of time in both of its numerator and denominator 
(respectively Q and a). As a result, the time derivatives of perturbations belong to the 
subleading order of the slow roll expansion, AQ = (A, ip, Sip) = 0(e) AQ. 

Thus, we arrive at the following list of correlators at late times. The potential type 
correlators read 

<^> = <*"> = "<^> = Sa^jfP^ 

^=2Sw f ' (8 ' 13) 



where the parameter / is given by eq. (|7.3| ), while the kinetic type correlators are negligibly 
small 

(AQAQ'} = 0(e) aH(AQAQ), 
(AQ'AQ'} = 0(e 2 ) (aH) 2 (AQAQ), 
( AQAQ" ) = 0(e) (aH) 2 ( AQAQ), 

( AQ'AQ" ) = 0(e 2 ) (aH) 3 ( AQAQ ) (8.14) 
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In view of these relations, the terms that give the leading contribution to radiation 
currents are exhausted by a small fraction of terms in the equations ( |6.27| )-( |5.28| ). They 
include only the potential type correlators and read 
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J if 



v 

x^ 3 



4(A 2 ) + 6(AV> - —(A8ip) + -(if>6<p) - 4 
2(A 2 ) + 6(AiP) - -(A6cp) + -{iP6cp) - ^(6cp 2 ) 



The resulting radiation currents, thus, equal 
Ay? 4 A 
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J (p 



4 8647r 2 |e| 2 ^' 
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2 864tt 2 |£| 



(8.15) 
(8.16) 

(8.17) 
(8.18) 



Similarly to the onset of inflation, eqs.( [7.8|) -( [7.9|) , they are strongly suppressed relative to 
the classical values, —V(<p) ~ — A<y9 4 /4 and —V^cp) ~ — A<^ 3 by the factor A/|£| 2 ~ 10 -10 . In 
absolute units, the energy density of the quantum mechanical mode is given by 



X 2 \l + 5\ mj 



54|^| 3 (16tt 



2^2 



(8.19) 



Note that, in contrast to the onset of inflation, this energy is negative. Apparently, this 
is a manifestation of the ghost nature of the invariant physical mode q, whose kinetic term 
enters the action (|5.18|) with the wrong sign. In the slow roll approximation with negligible 



ip a constant energy density ( J8.19 ) in view of eq.( |3.30| ) - the conservation law for radiation 
current - generates the effective equation of state of the inflaton excitation mode 



£ i + p i = 0, e q < 0, Ht > 1, 



.20) 



which in the absence of other sources would maintain the Anti-DeSitter spacetime. Similar 
equation of state at the onset of inflation corresponds to the DeSitter case, so that the 
inflaton quantum excitation undergoes a sort of phase transition reversing the sign of its 
energy density. 

In view of the relation between the radiation currents, j q N = ipj v /4: + 0(e) , and approxi- 
mately constant value of j q N) dj q N /dt = 0(e), the quantum rolling force in eq.( |6.35|) , 



F q 



1 dj q N 
2H dt 



0(e), 



(8.21) 



vanishes in the leading order of the slow roll expansion. Thus, similarly to the onset of 
inflation, at late times of inflation epoch F q does not qualitatively change the cosmological 
evolution. 



9. Conclusions 

We have developed a general framework for effective equations of inflationary dynamics 
in quantum cosmology and for their quantum Cauchy problem with no-boundary and tun- 
neling quantum states. This framework combines the Euclidean effective action method and 
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the method of direct quantum averaging for calculations of two distinctly different parts of 
radiation currents - contributions of the field theoretical and quantum mechanical (minisu- 
perspace) sectors of the system. We focus on the latter and show that its calculation is based 
on explicit physical reduction for the spatially homogeneous cosmological perturbation. Be- 
cause of the ghost nature of this perturbation, its effect is not related to the conventional 
analytic continuation from the Euclidean spacetime. Rather, in the model of strongly cou- 
pled non-minimal inflaton it originates from the quasi-gaussian state which incorporates the 
tree-level and one-loop effects on the DeSitter instanton. It is, thus, irrelevant to the De- 
Sitter invariant Euclidean vacuum and cannot be obtained by analytic continuation from 
the Euclidean section of spacetime. This means that the universality of analytic continua- 



tion methods of [44] should not be overestimated - they apply to spatially inhomogeneous, 
particle like excitation but may fail for minisuperspace cosmological modes. 

Unfortunately, the dynamical contribution of the quantum mechanical mode to effective 
equations turned out to be disappointingly small - it is strongly dominated by the effective 
rolling force (contributed on equal footing by the classical term and the one-loop term due to 
the inhomogeneous modes). The property of its strong suppression by powers of 1/|£| <C 1 
was actually conjectured in HJ, and now it is quantitatively confirmed. Thus, the inflaton 
mode cannot change the dynamical predictions in spatially closed model with strong non- 
minimal coupling. As a model of the low-energy quantum origin of the Universe only the 
tunneling state remains observationally justified, because the no-boundary wavefunction 
generates infinitely long inflationary stage. The role of this mode should not, however, be 
underestimated, because its effect is model dependent, and might be important in other 
models generating initial conditions for inflation [[5|. Moreover, the quantum inflaton mode 
simulates the DeSitter and Anti-DeSitter effective equations of state, e + p = 0, respectively 
at the onset of inflation and at late times. The sign of its energy density contribution can 
change depending on the balance of the potential and kinetic terms of this ghost mode. 
Therefore, it is not quite clear at the moment, what can the role of this mode be at post 
inflationary epoch. A natural question arises if this mode can be responsible for the present 



day observable acceleration of the Universe [46] as an alternative to quintessence [47] or be 



capable of inducing DeSitter-anti-DeSitter phase transitions in cosmology? This question is 



subject to further studies |45] 
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